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Abstract 

In stiff fiber-reinforced material, the high shear stress concentration occurs in the 
narrow region between fibers. With the addition of a small geometric change in cross- 
section, such as a thin fiber or a overhanging part of fiber, the concentration is sig- 
nificantly increased. This paper presents mathematical analysis to explain the rapidly 
increased growth of the stress by a small particle in cross-section. To do so, we consider 
two crucial cases where a thin fiber exists between a pair of fibers, and where one of 
two fibers has a protruding small lump in cross-section. For each case, the optimal 
lower and upper bounds on the stress associated with the geometrical factors of fibers 
is established to explain the strongly increased growth of the stress by a small particle. 

MSC-class: 35J25, 73C40 

1 Introduction 

In this paper, we concern ourselves with the high stress concentration occurring in the 
stiff fiber-reinforced composites when fibers are located closely. The primary investigation 
focuses on the case when a smaller fiber is located in-between area of two fibers, see Figure 
[T]and Figure [2] This paper reveals that, with the addition of a smaller fiber, the growth 
of stress is significantly increased: if the diameter d of the fiber in the middle is sufficiently 
small and the distance between adjoining fibers is e, then the stress blows up at the rate of 
in the narrow region, even though the blow-up rate has been known as as in the 
case of a pair of fibers. This means that the defect of fiber as a protrusion causes much 
lower strengths in composites than had been thought. To derive it, we estimate the optimal 
lower and upper bounds of the stress concentration in terms of the diameters of fibers and 
the distances between them. These bounds explain the dramatic change of the growth of 
stress when the diameter of the fiber placed in the middle is relatively smaller than other 
two fibers. 

In the anti-plane shear model, the stress tensor represents the electric field in the two 
dimensional space, where the out-of-plane elastic displacement satisfies a conductivity equa- 
tion, and the cross-section of stiff fibers corresponds to the embedded conductors. In this 
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respect, we consider the gradient of the solution to a conductivity problem to estimate the 
stress. Adjacent stiff fiber-reinforcement induces the high stress concentration in the narrow 
region between the fibers. This implies the blow-up of the gradient of the solution between 
adjoining conductors, seeOEHlEliaiinillllllS]. 

Meanwhile, the extreme conductivities are indispensable to the blow-up phenomena: 
when the inclusion's conductivity is away from zero and infinity, the boundedness of the 
stress function has been derived by Li and Vogelius [H] , see also [S] , and it was generalized 
to elliptic systems by Li and Nircnberg [TT] . In O H] , for conductivities including both 
bounded and extreme cases, Ammari et al. have established the optimal bounds of the 
gradient of solutions to the conductivity equation, when conductors are of circular shape 
in two dimensions, and the optimal bounds provides blow-up rate, where e is the 

distance between two conductors. Yun [T31[T3] has extended this blow-up result for the case 
of two adjacent perfect conductors of a sufficiently general shape in two dimensions. In Bao, 
Li and Yin's paper [7], it has been also investigated as the blow-up phenomena in higher 
dimensional spaces, also see [2 [T3]. They has also done a natural follow-up in [3 [S] that 
the blow-up rate known only for a pair of fibers is still valid for the multiple inclusions in 
any dimensions. 

In contrast, our paper witnesses an unexpected fact on multiple inclusions that the 
growth of stress can be significantly increased by a little geometric change of an inclusion, 
even though the blow-up rate is still e~^/^. 



For I = 1, . . . ,L, let Di be conducting inclusions in R^, that is cross-sections of stiff 
fibers. Then, under the action of the applied field H, the electric potential u satisfies the 
following conductivity equation: 




Figure 1: Case (A) and case (B) 




Figure 2: Case (C) and case (D) 



Au = 0, 
^ u(x)-ff(x)=0(|x|-i), 



u\dDi — C'l (constant), 



inM2\uf^,A, 
as |x| — !■ oo, 
for ? = 1, . . . , L, 
for / = 1, . . . , L, 



(1) 
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where H is an entire harmonic function H in and x = (.ti,X2). In this paper, we only 
consider the case of L = 2, 3, 

As we mentioned previously, for the two fibers with the circular cross-sectional shape, 
Ammari et al. [31 Hj obtained the optimal blow-up rate e^^^^ for Vm, and this result is 
extended by Yun [141 115j to general shaped fibers. Building on the prior results, we extend 
these into the the following interesting direction: we first consider the circular inclusions in 
Case (A) and Case (B), and second extend the result into general shaped ones in Case (C) 
and Case (D). In Case (A) and Case (B), we add a small circular inclusion between two 
others so that three disk centers are lined up in one straight line. The additional disk can 
be embedded disjointly from other disks, or it partially overlap one of two disks, and we 
formulate these two cases as follows. 

(A) One disk and a pair of partially overlapping disks, Figure [TJ there is a portion of disk 
protruding from one of circular inclusions, i.e., L = 2, and Di and D2 arc e-distanced 
domains defined as 

Di ^ BrAci) and D2 = B.Jcz) U 5,3(03), (2) 
where Bn (c; ) is the disk with the radius r.i and centered at c/ , and 

Cl = (-ri - |,0), C2 = (7-2 + |,0), and C3 = (rg + a+ |,0). (3) 

Here, 5^2(^2) is a small disk protruding from _B.r3(c3), and we assume 
5^2(02) n 5,-3 (C3) ^ 0, i.e., < a < 2r2, 
dist(Sri (ci), 5,3(03)) ~ r2 

and 

< e < 7-2 < 7'i ^ T-a- 

(B) Three disjoint disks. Figure [1] a small disk is disjointly embedded into the in-between 
area of two disks, i.e., L ~ and 

Di=Br,{ci), Z = 1,2,3, (4) 

where Ci = (-ri - C2 = (r2 + ^,0) and C3 = (r3 -I- r2 + ^ + 62,0). Hence, the 

distance between Di and D2 is ei, and the distance between D2 and I?3 is £2- Here, 
D2 is regarded as the cross-section of the thin fiber between a pair of fibers with the 
cross-section Di and £'3. Thus, we assume that 

< < r2 < ri ~ r3 for i = 1, 2. 



In both cases, the blow-up rate is remarkably increased due to the existence of -6,3(02) 
as follows: 

Theorem 1.1 (Case A: Protruding small disk) Let Di and D2 be defined as Then 
there is a positive constant C independent of e, ri, r2 and such that 
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>C- ^ 
ODi ri -I- rs Jr2 
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where u is the solution to ([T]) with H{xi,X2) ~ xi. As a result, by the Mean Value Theorem, 
there is a point xq in the narrow region between Di and D2 such that 



f 1 + rs Jrl ^ft 



For any entire harmonic function H , let u be the solution to (U]) with H . Then, there is 
a positive constant C independent of e, ri, r2 and r^ such that 

iv.(x)i < 

in the narrow region between Di and D2. 

Theorem 1.2 (Case B: Disjointly embedded small disk) Let Di, i ~ 1,2,3, be balls 
defined as Then there is a positive constant C independent of ti, €2, ri, r2 and rs such 
that 

-u >C ■ ^Vei, 

dD2 aoi ri + rs ^/r^ 

and 

rirs 1 

— u ^ u- 

dD3 

where u is the solution to ([T]) with H{xi,X2) = Xi. As a result, by the Mean Value Theorem, 
there exists points; xi in the narrow region between Di and D2; X2 in the narrow region 
between D2 and D3, which satisfy that 

\\7u(xA\>C^-^ fori^l, 2. 

n + rs Jr^ 



dD2 ri + rs Jr^ 



For any entire harmonic function H , let u be the solution to ^ with H . Then, there is 
a positive constant C independent of ei, €2, ri, r2 and r^ such that 

\Vu(^)\<C^^^^ for 1 = 1, 2. 

ri + ra Jr2 



in the narrow regions between Di and D2, and between D2 and D3, respectively. 

In this paper, we first estimate the lower-bounds in terms of the radii of inclusions. 
Based on this estimates we derive the remarkable blow-up rate increasing phenomena when 
a small conducting inclusion is located in-between region of two inclusions. This paper is 
organized as follows: In section 2, we explain the method to calculate the potential difference 
in the case of two disks. We then derive the lower bound of Case (A) in section 3; Case 
(B) in section 4. In the case of the upper bounds, the major part of derivation overlaps in 
Case (A) and Case (B). Thus, the derivation is presented in Subsection 14.51 Based on the 
similar derivation, we can also obtain the analogues of Theorem II . 1 1 and 1 1 . 21 for the inclusions 
associated by a sufficiently general class of shapes. 
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Analogues of Theorem 11.11 and 11.21 for a sufficiently general class of shapes 

The proofs of Theorem 11.11 and 11.21 arc flexible enough even though the results are restricted 
to circular inclusions. The estimates presented in Theorem 11.11 and 1 1 . 21 can be extended to 
the inclusions associated by a sufficiently general class of shapes. To consider a large class 
of shapes, we make the geometric assumptions more precise. To define Di, D2 and D3, 
we consider three domains ^'j-jgjit;' -^center ^^"^ -^left addition, we assume that 

V'right ■■ C\Si(0) K'\i?right' "^center : C\i?i(0) ^ R^XD^entev and y^j^ft ■ C\Si(0) ^ 
are conformal mappings in C^(C\i3i(0)) such that f'j.[gi^^{^) 7^ and f'^^^^iz) 7^ 

for z e dBi{0). Here, we do not distinguish from C. The regularity condition of these 
conformal mappings doses not allow non-smooth inclusions such as polygons, but Riemann 
mapping theorem yields a sufficiently general class of shapes: refer to Ahlfors [I]. Now, we 
consider the analogues of Theorem 1 1 . 1 1 and 1 1 . 21 for two cases as follows: 

(C) One domain and a pair of partially overlapping domains, similarly to Figure O there 
is a small portion of another domain protruding from a inclusion, i.e., L = 2, and Di 
and Z?2 are e-distanced domains defined as 

-Oi = ^left and D2 = (r2-Dcenter) U i'right' i^) 

where ''2^'center ''s times diminished domain of ^'center- suppose that D2 

is a connected domain, dist(£'i, 1)2) = dist(-Di, r2-DQQjj^tQj-), 

dist{Di,Dj.^^YiO - ^2, 
£>! C M- X M and Da C M+ x M. 
In addition, we also assume that 7'2 is small enough and 

< e < r2, 

and that the boundaries dDi, dD2 and dD^^^g^^ are strictly convex in the narrow 
region between Di and D2. 

(D) Three disjoint domains Di, D2 and D^, Figure O a small inclusion D2 is disjointly 
embedded into the in-between area of two other domains, i.e., L = 3, and 

-Oi = -^left' ^2 = 7-2£'center and = ^'right (6) 

where ''2 ^center ^^1*^ ''2 times diminished domain of ^center- We assume that 

Di and D2 are ei apart, D2 and D3 are £2 apart, and Di distances enough from D3 
that r2 is sufficiently small and 

< < r2 for i = 1,2, 

since D2 is regarded as the cross-section of the thin fiber, that the boundaries dDi, 
dD2 and dD^-^^,^ are strictly convex in the narrow region between Di and D2, 

C IR_ X M and U Z^a C E+ x M. 
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Theorem 1.3 (Case C: Protruding small lump) Let Di and D2 he defined as ([5]). Then 
there is a positive constant C independent of e and r2 such that 



dD2 



where u is the solution to ([T]) with H{xi,X2) = x\. As a result, by the Mean Value Theorem, 
there is a point Xq in the narrow region between Di and D2 such that 



1 1 



|Vw(xo)| > C 



For any entire harmonic function H , let u he the solution to ([T]) with H . Then, there is 
a positive constant C independent of r2 and e such that 



|Vu(x)| < C 



1 1 



in the narrow region between Di and D2. 



Theorem 1.4 (Case D: Disjointly embedded small inclusion) Let Di, i ~ 1,2,3, be 
balls defined as Then there is a positive constant C independent of r2, Ci and 62 such 

that 

dD2 dDi Jrl 



dD2 Jr2 



where u is the solution to ([T]) with H{xi,X2) = xi. As a result, by the Mean Value Theorem, 
there exists points; xi in the narrow region between Di and D2; X2 in the narrow region 
between D2 and D3, which satisfy that 



|Vm(x,)| >c 



1 1 



r2 Ve -t 



for i = l, 2. 



For any entire harmonic function H , let u be the solution to ([T]) with H . Then, there is 
a positive constant C independent of r2, ei and £2 such that 



|Vw(x,)| < C 



1 1 



for i^l, 2. 



in the narrow regions between Di and D2, and D2 and D3, respectively. 

Wc derive the lower bound of Case (C) in section 3; Case (D) in section 4. In the case of 
the upper bounds, the major part of derivation overlaps in Case (A), Case (B), Case (C) 
and Case (D). Thus, the main idea is presented in Subsection [ 
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2 Preliminary 



2.1 Calculation of the potential difference 

We explain the main idea to calculate the difference of potential between two adjacent, 
possibly disconnected, conductors. 

In this section, differently from ([T]), Di, i = 1,2, could be also the union of two disjoint 
domains. Define u as the solution to ([T]), where it is assigned one constant value throughout 
Di even when Di is disconnected. Now, define h as the solution to 



Ah^O, iiiW\{DiU D2), 

h = 0(|xp^), as |x| 00, 

h\dDi ~ ki (constant), for i = 1,2, 

/ d,,hdS^{-iy, fori = 1,2, 



(7) 



where is the outward unit normal vector of R"\(Z)i U D2), i.e., directed inward of Di. To 
indicate the dependence of u and h on Di and D2, we denote them as 

w = $[Di,D2], (8) 

h^^[Di,D2]. (9) 
The potential difference of u in Di and D2 is represented in terms of h as follows. 
Lemma 2.1 (fW) 



dD2 



Hd^h dS - 



Hd^h dS. 



dD2 



(10) 



The lemma above can be derived by the Divergence Theorem, see 



2.2 Two disks in 

Using Lemma I2.H we can easily calculate the potential difference m|d2 ^"l-Di of the solution 
u to ([T|) when 

=B,,(ci) andZ?2 = S,,(C2), (11) 

where Ci = (— ri ~ |, 0) and C2 ~ {r2 + f , 0). 

Let Ri be the reflection with respect to Di, in other words. 



i?,(x) 



^(x-c,) 



Ci, i = 1,2, 



and Pi G Di be the fixed point of Ri o i?2, then i?2(Pi)(=: P2) is the fixed point of R2 o Ri, 
and 



Pi 



V2 



rir2 



n + r2 

Moreover, we can easily show that 



V^+O(e),0) and P2 =(\/2 



rir2 

ri + r2 



V~e + O{e),0 



^[D,,D2] = — (log|x-p,|-log|x-p2|). 



(12) 
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By an elementary calculation, it can be shown that the middle point ^'-^^^ exists between 
two approaching points (— |, O) and (f , O) . Applying the middle point property to estimate 
for D2] (±|, 0), we can get the following lemma. 

Lemma 2.2 There is a constant C > independent of e, ri and r2 such that 



for small e > 0. 

From Lemma we calculate the potential difference of u. 

Lemma 2.3 Let H{xi, X2) he an entire harmonic function. The solution u to (QJ) where 
L ~ 2 and Di, I = 1,2, are given as (jlip satisfies 

(13) 



2V29.,i/(0,0)y^V^+O(e). 



Remark 2.4 Referring to the mean value theorem, there exists a point X2 between dDi and 
dD2 such that 

\Vu{x2)\ > 2x/2|9,,ff(0,0)L/^^^. (14) 

V ''1 + ''2 

for any sufficiently small e > 0. Moreover, as a result in J4^, there is a constant C indepen- 
dent of e, ri and 7^2 such that 



I 7*1 r2 1 



+ r2 i/e 
where fl = i?4(ri+r2)(0, 0). 

3 One disk and a pair of partially overlapping disks 

In this section, we consider two e-distanced domains Di and D2, see Case (A) at Figure [1] 
where 

L>i = Bri(ci) and D2 = Br,{c2) U ^^3(03), 

where 

ci = (-ri - -,0), C2 = (ra + -,0), and cg = (rg + a + -,0). 
Here, -8^2(^2) is a small lump of _Br3(c3), and we assume 

^^2(02) n 5^3(03) ^ 0, i.e., < a < 2r2, 

and 

and < e <C r2 <C min(ri, r^). 

Define 

h ^ '^[Di, D2], and u = <^[Di, D2], (15) 

and 

h, = ^[Di, Br^ (c,)], and u, - $[^1, S,^ (c,)], j = 2, 3, (16) 
where ^ and <f> defined in section [TT] 
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3.1 Properties of h and hj, j = 2,3 
Lemma 3.1 Let h = '^[Di, D2], then we have 

dJi{x) = 0{V^), xedD2\Br,{c2). (17) 

Proof. Set 

B, = BrXci), fori = 1,2, 3. 
We choose a smooth domain 51 as follows: 

r2 C (B2 U B3), B-idfl 
dh\B2 = dBs \ B2 

{dnndB2)\dB3 = (^6,0) 

Let h = 'S[Bi,n]. Then, we consider V defined in \ (Bj U B2 U B3) as follows: 

h\dBi -h\g^ 



Then, it follows that 

V\dBi = y\on\B2 = ^ constant. 

Since /i|a(B2uB3) > ^lasi £^nd /ilg^ ^Is-Bu the minimum ofV attains on dQ\B2 and dBi. 
Thus, we have 

d,h - ^'^f ^ ~ ^l^i^-^^^) d~h < on U (dn \ B2). (18) 

h\dB, - /llgn 

By the integration on dBi, we have 
Using the bound (fTH|) once more, we have 

< on af2 \ B2 = \ B2. 

The domain 17 is smooth so that we can use the method presented by Yun [TU [TS] . Then, 
up to a conformal mapping to a circle, dyh is bounded by constant times the Poisson Kernel 
with respect to a interior point ^/e distanced from the boundary (refer to the inequality (9) 
in [H]). Note that 9.83 \ B2 distances enough from (e, 0). Thus, we have 

duh< dJi\Q^^ < CV^ on dB^ \ B2. 

Therefore, we have completed the proof of the lemma. □ 
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d^x) < Md^hsix), X e dDi 



M 



Lemma 3.2 

where 

Proof. Define 

W^h-Mhs, ii-iR^ \iDiU D2). 
Since h is constant on dD2, M > 0, and ^,3 takes it's maximum on 95^3 (ca), 



h 




dD 


1 




93.3(03) 




dDi 



(19) 
(20) 



and 



w 


-w 




= ~1 




dD2\Br^(c2) 


dD2\B,3 (03) 


w 


- W 


= h 


- h 




9-D2\Br2(c2) 


dDi 


dD2 



dDi 



?Br-3(C3) 
Af(/l3 



dD2\B,^ (C3) 



) <o, 



dDi 



)=0. 



Therefore, W takes its minimum on and 

dyW < 0, on dDi. 



Lemma 3.3 



Proof. Note that 



- h 


= h2 


dD2 


dDi 



dBr2(c2) 



dDi 



0{e). 



□ 
(21) 



and 



/ d,ih - h2) dS = 0, 

JdDi 

I d^{h~h2)dS^f d^^hdS- [ d^h2 dS - [ 

JdD2 JdD2 JdB„(c2) Jd 



duh2 dS 



d(D2\B^^{c2)) 



1-1-0 = 0. 



With the fact that /ilauj and h\go2 arc constants and the (exterior) Divergence Theorem, 
we have that 







Hence, 



dDi 



dD 



[ d^{h-h2)hdS+ [ d^{h-h2)hdS 

JdDi JdD2 



{h - h2)d^h dS + / {h-h2)d^hdS. 

dDi JdD2 



hd^h dS + / hd^h dS 

2 JdDi JdD2 



h2d^hdS+ / h2d^hdS 

ODi JdD2 



dDi 



+ / ih2~h2 

aS,-2(c2) JdD2 



)d^h dS 
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From ([T^ . there is a constant C dependent of a, see ([3]), such that 



(/i2-/i2 ){x) <C^fe, forallxe9D2\S,.,(c2). 



as^2 ('=2) 

Therefore, with ([TT]) as well, we obtain pT|) 
3.2 Proof Theorem [TT] 



□ 



Let H{xi, X2) = xi and be the unit normal vector of M^\(Z)i U D2); i.e., directed inward 
to Dj, i = 1, 2. Remind that we defined 



h = ^Di,D2l andu = <i>[Di,D2], 



and 



(22) 
(23) 



hj = -^[Di^Br^ {cj)l and Uj = S,.^ (c^)], j = 2, 3, 

where ^' and $ defined in section [TT] 

Note that < 0, iJ < on dDi and H > on 91)2, and, as a result, from Lemma 

dD2 

12.11 we have 



dD2 



= / {d^h)HdS+ j {d^h)H dS 



> / Hd^hdS. 

JdDi 



Applying the lemma [221 Lemma 



becomes 



dD2 



> 



> 









1 




05^3(03) 




dDi 


h2 


as,.2(c2) 


-h2 


dDi 






(C3) ~ 


h3\ 



Hd^hs dS 



0{e) 



V2 



\ 



It follows from Lemma [2?2] that 



as,2(c2) 



dDi 



> C 



ri + r2 
rir2 



and 



Therefore, 



3.3(^=3) 



dD2 



<cJ^l±^V^ + 0(r2). 
dDi V nrs 



dDi n + rs Jr2 



(24) 



This proves Theorem 1 1.1 



□ 
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3.3 Proof Theorem Q 



Wc consider the general shaped domain in Theorem 11.31 But. we take an advantage of 
the properties of circular inclusions. To make a connection between circular domains and 
general shaped domains, we need to establish the monotonic property of 5* as follows: 

Lemma 3.4 [Monotonic property of 5*/ Let Da, Db, Da and Db be domains. Assume 
that ^ ^ 

Da C Da and Db C Db- 

Then, we have 



0<^DA,DB]\gS^-MDA,DB]\gg^<nDA,DB] 



ODe 



^Da^DbW 



dDA- 



Proof. Without any loss of generality, we consider only the case of Da = Da- Let 

G^-^[Da,Db]- AmDA^De] 

where 

M 



^[Da^Db] 


^^^-^Da.Db] 


dDA 


^Da^Db] 


^^^-^Da^Db] 


ODa 



The minimum of G attains on ODa- By the Hopf's Lemma, we have 

a^G < on ODa. 
Integrating d„G on ODa, we have — 1 + 7\/ < 0. Therefore, we have 

Q<nDA,DB]\Q^^-nDA,DB]\Q^^<nDA,DB]\,^^-nDA,DB] 
Repeating the same argument again, we can obtain the disable inequality. 



□ 



Applying -Djg^^, 7'2^'center ^^"^ ^right instead of Bri{ci), i =■ 1,2,3, to the argument 
presented in the proof of Theorem 11.11 we can obtain 



^'■^^center 



dD 



> C- 



(''^ ^center) 



\dD 



left 



0(6) 



left 



^^right 



h3 



OD 



left 



when H{xi,X2) = xi. Here, hi = ^'[-Dieft> ''2^'center] and /12 = *[-Dleft' ^rightl- 
It follows that from Yun [HI [H] that 



- h. 



dD, 



right 



dD 



left 



To estimate 



91-2 -D 



center 



dD 



left 



, we choose two disks Sj^f^ and B(,q^iq^ containing 



^left a^*^ ^center such that the distance between i?igf^ and f2^center e. Using Lemma 
Ell and [13 we have 

^2 — ^2 ^ 

^('■^■"center) ^-"left 
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Note that Di = D^q^-^, D2 ~ '"2 -D^gj^^^gj. and = -D^jgi-^^^ in this theorem. Therefore, 



dD2 



dDi Jrl 



This proves Theorem 11.3 



□ 



4 Three disjoint smooth domains 

We consider three disjoint inchision case, see Figure [T] and [3 a small one is disjointly 
embedded into the in-between area of two others, and prove Theorem 11.21 and 11.41 We 
assume that Dx and Di are closely spaced with the distance ei, and and -D3 are closely 
space with €2, but D\ and D3 are not close, and that Di, Di and Dj, have the boundary 
regularity given in Theorcm ll.4l 



4.1 Solution representation of u 

Let H'^ be a harmonic function outside of vS^^^Di and have the same constant value in 
Uf^jDi satisfying that 



Ai?'^ = 0, in K2\u3^iA, 

H^(x) - i/(x) = Cl(|x|-^), as|x|^oo, 
i7^|y3_^g£,. = Cu (constant). 



(25) 



Since H'^ — H is harmonic at infinity, H'^ — H attains maximum only at the boundary points 
of -Di, i = 1, 2, 3. To make H'^ — H attains zero at infinity, Ch should satisfy 



II^IIl~(u3^iD,) - ^-f^ -l|-^IL~(uf^iD,)' 

Moreover, iJ^ satisfies Y^i=i Ian ^j^^^ d,S = 0. 
The solution u to ^ is represented as 

m(x) = H^ix) + Ci/li(x) + C2/l2(x), 

where 



(26) 



(27) 



and 



ci 

C2 



hi = *[A,(^2Ui?3)], h2^'<S/[{DiUD2),D3 

/ 1 I a u \ / 



-1 



V 



dDi 



J 



d.H' dS 



\ 



dDi 



V 



d.H' dS 



(28) 



dD2 



J 



dM dS / d^h2 dS 

ldD2 JdD2 

where ^ is defined as ([7]) and The equality (pS)) is from the integration of d^u on dDi 
and dD2. 

Applying the upper bound on the gradient of solution without the potential difference 
among the boundaries to conductivity equation derived in Bao et al. [7] , we can show that 
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VH'^ does not blow-up (also refer to [II]). Using Lemma [4.31 in the following section, we 
have duhi dS = 1 + 0{y/ei). This implies that 



ci 

C2 



-1 
1 -1 



dDi 



\ J 00-2 



d^H"' dS 
d^H" dS 



\ 



Thus, the coefficient i = 1,2, is bounded independently of ei and £2. Therefore, the 
blow-up rate of Vit essentially relies on V/i^. In this respect, we consider the properties of 
hi in the following section. 

4.2 Properties of hi and h2 

We build the optimal bounds of u based on (|27p ; it is essential to drive properties of hi and 
/i2 in the narrow regions between inclusions. Let hi and /12 be as follows: 

hi = «'[L'i,(7:>2U A3)] and /12 = [pi U L>2), £'3] • 
Proposition 4.1 There are the following estimates for hi and /i2. 

(i) In the narrow region between Di and D2, we have 

Whi = O(^) and V/i2 = 0{y/^). 



(ii) In the narrow region between D2 and D3, we have 



Vhi = O(Vei) and V/i2 = O 



(Hi) 



and 



hl\dD2\JdD3 — hi\dDi — 

/«2|a_D3 — /j2|a_Diuar>2 — \/^- 



Proof. We consider V/ii. By Lemma l4~2l and [4~4l we have 

> d^hi > Ca^*[I?i, Di] on dDi 



and 



and by Lemma l4.3 



< d^hi < d^'^[Di,D2] on dD2, 



0< \d,hi\ <C^i on 
where Di is defined in Lemma Without any loss of generality, we assume that 



(-^,0) ea/^i, (^,0) ^dD2 anddist(Z?i,D2) =£1. 
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Let 

p(x) = log |x - (y?r, 0)1 - log |x + (ViT, 0)|. 
Referring to the inequality (9) in [T5|, there is a constant Ci such that 

< I V/ii I < Ci |Vp| on d{Di U L»2 U D3). 

Regarding {xi,X2) as a complex number z = .ti + 2:21, we consider 

dihi{z) - d2hi{z)i 



p{z) 



Ci{dip{z) - d2p{z)i)' 



Then, p{z) can be extended to 00 as an analytic function. From definition, \p{z)\ < 1 on 
dDi U 9Z?2 U dD^. By the maximum principle. 



|p(z)| < 1 mC\{DiUD2UD3). 



\\/h\ < Ci\Vp\ in R2 \ (^Di UD2UD3 



Thus, we have 



Therefore, V/ii = ^(^^j in the narrow region between Di and D2, and V/ii = 0{^/ei) 

in the narrow region between D2 and D3. Similarly, we have V/7.2 = '^(^"^^) ^^'^ narrow 
region between D2 and D^^ and V/12 = O(y^) in the narrow region between Di and D2. 
We have proven (i) and (ii). 

The estimate (iii) is presented by Lemma 14.41 

□ 



Lemma 4.2 We have the following properties: 

< hi\9D2 - hi\9D, <*[i?l,£'2] 



(^) 



dD2 



nDi,D2] 



aDi 



< d^hi < 0^-^101,02] on dD2. 



Proof. Let 



and 



M = 



^[Oi,{02UD3)] 


(i^2Ui?3)] 

dDi 


9D2UOD3 




^■^[0^,02] 

dDi 


dD2 



G(x) = ^Oi, {D2 U DsM^) - {O2 U O3] 

-M('f[Oi,02]{x)-nOi,02] 



dDiUdDa 



dD2 



Then, G = on dOi U d02, and G > on dO^. By Hopf's lemma, 

9^G<0 onaz?!. 
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This means that 

d^hi < Md^^[Di,D2] on dDi. (29) 

Note that hi = ^'[£'1, {D2 U D^)]. By integrating G on dDi, we have the inequahty (i). 
On the other hand, by Hopf 's lemma, 

< on dD2. 

This means that 

d^hi < Md^'9[Di,D2] on dD2. 
From the inequahty (i). A/ < 1. Therefore, we have (ii). 



□ 



Lemma 4.3 There is a constant C such that 

< d^hi < Cy/^i on dD^. 

Proof. We use the method similar to Lemma IT^ Let 



^[Di,{D2UD3 



M = 



dDi 



■^[Du{D2\JDi)] 



dD2UdD3 



dDi 



and 



G(x) = {D2 U I?3)](x) - nOi, {D2 U D3] 

~m(^[Di,D3]{x)-^Di,D3] 



dDiUdDa 



dDs 



Then, G = on dDi U dDs, and G > on dD2. By Hopf's lemma, 

d^G < on dDs. 

Since hi = ^[Di, {D2 U D3)], this inequality means that 

< d^hi < Md^'i'iDi^Ds] on dDs. 

Now, we estimate the gradient of M'^[Di, D^]. To do so, we consider the potential 
difference between dDi and dDs as follows: 



M^[Di,Ds] 



M%Di,Ds 



dDi 



h 



h\dDi 



h\dD2 - h\oDi 



dD2 



^[Di,D2] 



dDi 



<C^i 

The last inequality above was proven by Yun in his paper [141 115] . since 5'[Z?i,Z)2] is only 
for two domains. Note that Ds is not close to D2. Owing to the method in Bao et al. [7], 
we have 

\\d,AmDi,Ds]\\L^^an,)<C^i. 
Therefore, we can obtain the result. □ 
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Lemma 4.4 Let D4 is a disk containing D2 and D3 with 

dist{Di,Di) = dist{Di,D2). 
(i) There is a positive constant C such that 

> djii > Cd^^[Di,Di] on dDi. 



(Hi) 



hi\dD-iyjdD3 — hi\QDi — Vei- 



dDi 



Proof. To prove (i) and (ii), we use the same derivation to Lemma [4?2l So, we set 



M = 



dDi 



dD2UdD3 



dDi 



and 



G(x) = ^[Di, {D2 U D^)]{^) - {D2 U D^)] 



dDiUODa 



dDi 



Then G 



ODi 



and G < on dD^. By Hopf's lemma, we have 

> on dDi. 



By the integration on 9Z?i, we have (ii) and M < 1. Therefore, the inequahty cJ^G > can 
also yield (i). 

From (i) of Lemma 14.21 and (ii) in this lemma, we have 



and 



hi\oD,-hi\9D,>MDi,Di] -^Di,Di] 

dDi 



hl\9D,-hi\gD, <MDl,D2] -^[Di,D2] 

dD-2 



dDi 



dDi 



The potential '^[Di^Di] (i = 1,4) is only for two domains and thus, its difference between 
Di and Di {i ~ 1, 2) was already estimated in Yun [TH[TS] as follows: for i = 1,2, 

^[Di,D,] A] ^V^. 

dDi ODi 



Therefore, we have (iii). 



□ 
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Lemma 4.5 We have 



Hd^hi dS 



Proof. Without loss of generality, we assume that 

[-J^O) e dDi, (^,0) G dD2, distpi,Z?2) = ei and (-1,0) e D2. 



We consider H as follows: 



X2 



x-(l,0)|2- 



It follows from the Divergence Theorem that 

X2 



dDiUdDi 



uaD3 |x - (1,0)1- 



:d^hdS = 



dDiUdD2UaDa 



X2 



''V|x-(1,0)P 



hds = 0, 



|x-(l,0)P 



0(|x| ^) as |a;| 00. Hence, we have 



Hd^^hi dS 



dDiUdD2 



Hd^hi dS + / Hd^hi dS. 



We first consider jgj^i^yjQjj^ Hdyhi dS. By Lemma W?^ and WM, we have 

> djii > Cd^'^lDi^Di] on dDi 

and 

< d^hi < a^*[Di,D2] on dD2. 
From definition, d2H = 0. Hence, we can use Lemma 3.2 in [T^ so that 



Hd^hi dS 



dDi 



and 



Hd^hi dS 



< C 



< 



H^[Di,Di] dS<Cy/e{ 



dD2 



dD2 



H'9[Di,D2] dS<Cy/^. 



We second consider Hdi,hi dS. By Lemma [4.31 '^e can have 



Hd^hi dS 



Therefore, we have done it. 



□ 



Remark 4.6 We draw attention of readers to the independent work of Bao, Li and Yin in 
fSj/ and JSjj. Bao et al. have shown that the blow-up rate know only for a pair of inclusion is 
still valid to the multiple inclusions cases. As a byproduct of our work, the blow-up rate of 
the gradient for three inclusions is established in Theorem \4. 7[ 
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Theorem 4.7 Let Di, D2 and be as assumed in the beginning of Section^ Note that 
D2 is not assumed to be smaller than the others. 

(i) Optimal upper bounds: For any entire harmonic function H[xi,X2), we have the fol- 
lowing: in the narrow region between Di U D2, 

\Vu\<C-^, 



and, in the narrow region between D2 U D^, 

'Vitl < C- 

£2 



(ii) Existence of blow-up: Without loss of generality, we assume that 

(-|-,0) <EdDi, (|-,o) <=,dD2 and dist{Di,D2)-^ei. 
For H{xi,X2) = Xi, there exist Xq in the narrow region between Di and D2 such that 



and, similarly, there is a linear function H{xi,X2) with yo between D2 and D3 such 
that 

\yu{y,)\>C^. 



Proof. From Subsection 14.11 wc have a representation (pS)) for u and the coefficient Ci , 
i = 1,2, is bounded independently of ei and £2. Proposition 14.11 yields the upper bound of 
Theorem HTl 

Now, we consider the existence of the blow-up. Using the result of Subsection 14. II again, 
we have a constant C independent of e such that 

ll"llL'°(u-^^iaD,) < C'!l^llL°°(uf^ini)- 
Applying the Green's identity to J^js q^, udjii dS, we have 

Hd^hi dS= [ ud^hi dS 

~u\ 



^ ^ JaD2 ^ JdDa 

1- [ d^hi dS)+u\ ( / d^hi dS). (30) 



By Lemma 14.31 we have 

n , — III 

\dD2 \dDi 



wU„ -uU„ < C^, (31) 
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where the constant C above depends on ||-f^||L=o(y3_^3^.-). Similarly, we have 



(32) 



Hd^hi dS 



Using ([50]) again, we have 



> / Hd,,hi dS. 

JdDi 

The last inequality can be derived from the fact that H > Q on dD2 U dD^. 
To get the last inequality above, we took an advantage oi H = xi. By 



(33) 



hi 



dDi 



hi 



dDzUdDs 



dDi 



d^'9[Di,D2] < on dDi. 



dD2 



By (iii) in Proposition 14.11 we have 



hi 



dDi 



hi 



and 



^[Di,D2] -*[i?l,i?2] 



O-D2 



The inequality ([55]) implies 



9152 



> 



By the Mean Value Theorem, we have the desirable lower bound in the narrow region 
between Di and D2. Similarly, we can also obtain the other lower bound. □ 

4.3 Proof of Theorem D 

We derive the optimal bounds of the gradient of the solution to ([T]), when there arc adjacent 
three disks: 

Di^Br.ici), Z-1,2,3, (34) 

where Ci = (— ri — ^, 0), C2 = (r2 + ^, 0) and C3 = (r's + ^2 + y + ^2, 0). As defined before, 
hi = ^Di,{D2UD3)]. Let wi ^ ^Di,D2]. 
We begin the proof by showing that 



By the monotonic property of Lemma 13.41 we have 



Considering 



hi 



dD2 



Wi 



dDi' 



(35) 



hi 


dD2 


~hi 


dDi 


Wl 


dD2 


- Wl 


dDi 



Wl, 
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we can obtain, from the Hopf 's Lemma, 



hi\ 


9-D2 


-hi\ 






\dD2 


- Wi 


\dDi 



dMS < -^{^ / d.widS. 



By Lemma l4.3[ we have 
Since /a£),uar>3 d^hidS = 1, we have 



/ dMS - O(V^). 

JdD:i 



dDi ■ 



Therefore, we can obtain (1351) . Owins; to the estimate for wi L„ — wiLr, in Lemma 

^ ^ ^ \a-U2 \o-Ui 

we have 



^AaD2 - ^AoD, - V ^7^^- 



Let W2 — ^'[DiiDs]. Considering 

hi - I i i I W2, 

from the Hopf 's Lemma, we obtain 

duhi < ( '.^^ I d^W2 < on az?!. 



hi\ 




-/ill 


dDi 


^2 


a£'3 


— W2 


\dDi 



hi\ 


IdDs 




dDi 


W2\ 


IdDs 


- W2 


\dDi 



Here, we estimate the coefficient in the right hand side. Note that h-\\„^ = h-\\„^ . Thus, 
we have 

Since r2 ^ ri and r2 ^ J's, we also have 



This imphcs that 



Therefore, we have 



' 1-1+7-2 

< " ^ 9^W2 < on dDi. 

rirs V ^ 



/" HdMS > J '"^ ^ '"^ — / Hd^W2dS 
JdDi \l n+rs r2 JoD, 

\l ri+rs r2 \ ri+ ra 

>^-^^V^i>o. 37 
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Owing to (1201), (EB) and ([32]), we have 

[ ud^hi dS = (l- [ d,M ds\ fuL 

Juf_,dD, \ JdD3 J \ ^ 



IdDi 



dyhi dS j (u 



'dDs / 



Therefore, we have 



\aD2 "iar> 



\ [ Hd.hi dS + O(y^) (0(^^^) + 0{^)) 

^ JdDi 



2 

> 



+ ^'3 V^2 



Therefore, we have completed the proof. □ 
4.4 Proof of Theorem 11.41 

We pursuit the proof of Theorem 11.21 taking an advantage of the monotonic property of 
Lemma [3.41 The domains Di, D2 and D3 are as assumed in Theorem 11.41 As assumed 
before, hi ^ '^[Di, {D2 U D3)]. Let wi = ^'[L'l, £'2]. By the same way as Theorem [TT^ we 
have 

"^ilsD, -"'ilaDi -'^ilai?. -'^ilaiJi- 



Here, we use the monotonic property of Lemma 13.41 to estimate the difference between 
domains. Choosing two pairs of proper disks containing Di and D2, and contained Di and 
D2, respectively, we can obtain 

under the assumption that r2 is small. 

Let W2 = ^'[Di,£)3]. Choosing two pairs of proper disks containing Di and D^, and 
contained Di and D3, respectively. Then, we have 



W2\g^^-W2\g^^^^r2. 



By the same argument as Theorem II. 2 [ we have 



SDi V 
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Note that Di C ]R_ x M and D2UD3 cR+ x R. Continuing to follow the proof of Theorem 
11.21 we can obtain 



\dD2 



Therefore, we have done the proof. 

4.5 Derivation for the optimal upper bounds 

We consider the optimal upper bounds presented in Theorem 11.11 11.21 11.31 and 11.41 These 
proofs have essential thing in common. In this respect, we prove only the optimal upper 
bound presented in Theorem 1 1.21 As have assumed them before, we set 

hi = ^[Di,D2U D3] 
h2 = -^[01,02] 

/14 = ^[Di,D4]. 



Here, the domain D4 is given in Lemma |4.4| which is a disk containing D2 and with 

distpi, 1)4) = dist(i:'i, Liz), 

and the diameter of 1)4 is in proportion as r^, because r2 is sufficiently small. Then, we 
compare hi with h2, /13 and /14. The proof of Lemma 14.21 contains 



fhi 



< d^hi < 
the proof of Lemma 14.31 yields 
< d^hi < 



d{D2UD3} 



hi 



dDi 



h2 


-h2 






dD2 


dDi 



dyh2 on 9D2, 



(hi 



a{D2UD3) 



hi 



dDi 



V 



h3 


-h3 








dDi 



dyh-i on 9D3 



and the proof of Lemma 14.41 implies 

(hi 



< -d^hi < - 



a{D2UDs) 



hi 



V 



/14 


— /14 











dyhz on dDz- 



In the same way as Lemma [4.51 we can consider H by choosing the point in D^. In this 
respect, without any loss of generality, we can assume that 

9,,il(0,0) = 0. 

The reason why we assumed above is because the integration representation for the potential 
difference is not good enough, refer to [13]. The geometrical assumption of Case (B) implies 
that Di and D2 U -D3 are separated by Xi = and they are approaching to (0, 0). 
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Therefore, by the proof of Theorem II . 2 1 and Lemma [^751 we have 



\dD2 



IdDi 



+ o{^) (O(v^) + o(v^)) 

Hd^hidS 
hi 



'hi 



< 



d{D2UD3) 



dDi 



h2 


^h2 






dD2 


dDi 



hi 



d{D2UD3) 



hi 



dDi 



\ 


h3 


-h3 


dDi 


(hi 


d(D2UD3) 


hi 


dDi 


\ 


hi 


— hi 


dDi 



/_nT2_ 

T"! + f 2 ^ 



'— — r2 

n + rs 



ri + ^3 



-ei 



and 



hi 



-hi 




-h2 




a(r>2UAi) 




dD2 


dDi 



Here, note that the radius of D4 can be choosen between |r3 and 2r3. Lemma 
that 



imphes 



IdDi 



< 



ri + rs Jri 



\dD2 



\dDx 



Therefore, we estabhsh the optimal upper bound for 

Based on this, the optimal upper bound on the gradient of u in the narrow region be 
obtained. Here, the main idea to get the gradient estimate from the potential difference has 
already been presented by Bao et al. (Theorem 1.3, Lemma 2.2 and 2.3 in 0), and has been 
modified to fit our problem by Lim and Yun in |13) . Thurs, we give a brief description on 
the method. We choose a large domain Dq containing Z?i, Z?2 and D3, where 9Z?o is at a 
sufficient distance from Di, D2 and D3. Then, u can be decomposed as follows: 

u = Co + Wo + Civi + C3W3 

where for 1 = 0, 1,3, Vi is a harmonic function in Do \ (-Di U D2 U £'3) with the boundary 
data 

Vi — Soj on dDj for i = 1, 3 

and 

'Vq = ^ijU on dDj 
for any j = 0, 1, 2, 3. Thus, the constants Ci and C3 keep 

rirs 1 



\Ci\< 



and 



IC3I < 



n + rs y/r^ 

rir^ 1 
ri + rs Jt^ 
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To estimate Vwq, we consider a harmonic function p in Do \ {Di U D2 U D3) with the 
boundary data 

p = Soj on dDj 

for any j = 0, 1, 2, 3. By comparing with the harmonic function pi in Dq \ Di with pi = 
on and pi = 1 on ODq, the Hopf's Lemma yields 



I1Vp||l=°(_Do\(OiUD2UD3)) 

< max{||Vpi||Loc(a_Di), ||Vp2||L-(aD2), II Vp3||l-(sa3)> II Vp||l-(s_Do)} < C*- 

Applying the Hopf's Lemma again, we can have that the gradient of vo is bounded indepen- 
dent of ei, refer to Lemma 2.2 in |7j. 

We estimate CiVwi in the narrow region between Di and D2. Since vi is constat on the 
boundaries and the boundaries is smooth enough in the narrow region, the proof of Lemma 
4.3 implies that vi can be extend into the interior areas of Di and D2 by the distance almost 
e from the boundaries in the narrow region, independently of ri and r2. By the gradient 
estimate for harmonic functions allows 



in the narrow region between D\ and D2. Note that the inequality above is a local property 
independent of choosing Dq ■ 

Now, we consider CaVi'a in the narrow region between Di and D2. Let p be a harmonic 
function in Dq \ {D2 U D3) with the boundary data 

p = on d{DQ U D2) and p ~ 1 on dD^. 

By the maximum principle, we have 

0<V3<p in Da \ {Di UD2U D3) 

Considering the standard estimate for ^'[-D2, D3], we can obtain 

\C3V3\<C^. 

Similarly to the estimate for CiVwi, the gradient estimate for harmonic functions yields 



in the narrow region between Di and D2. 

Therefore, we can obtain the desirable upper bound. Here, it is noteworthy that the 
upper bound is dominated only by the estimate for C2Vi'i, which is independent of choosing 
Dq. In this respect, the constant C of the upper bound in Theorem 11.21 is independent of 
^'i, r2, rs, £1 and £2- □ 
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